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Abstract. We investigate the possible existence of non-topological solitons in string- 
like theories, or in other completions of Einstein theory, by examining a simple extension 
of standard theory that describes a non-linear scalar field interacting with the Einstein, 
Maxwell and Weyl (dilaton) fields. The Einstein and Maxwell couplings are standard while 
the dilatonic coupling is taken to agree with string models. The non-linear scalar potential 
is quite general. It is found to be impossible to satisfy the dilatonic boundary conditions. 
Excluding the dilaton field we find a variety of solitonic structures differing in ways that 
depend on the non-linear potential. In general the excited states exhibit a discrete mass 
spectrum. At large distances the gravitational field approaches the Reissner-Nordstrom 
solution. 
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I. INTRODUCTION 



Since the "black hole solitons" of string and supergravity theories carry horizons and 
central singularities, they differ from regular solitons which are singularity-free. It is natu- 
ral to ask if these generally covariant theories also admit regular non-topological solitons, 
similar to the field structures that exist at the level of special relativity, 1 ' 2 since it is known 
that non-topological solitons also exist at the general relativistic level. 3 ' 4 ' 5 It is possible, 
however, that there are generic features of supergravity and of other fundamental com- 
pletions of the Einstein theory that preclude the existence of everywhere regular solitons 
that are non-topological. In particular, this role may be played by the dilaton field, which 
arises as a consequence of Weyl rescaling, couples to the antisymmetric tensor fields, and 
replaces the numerical coupling constant. 

As these theories are derived by compactification from higher dimensional formula- 
tions, they characteristically contain multiplets of nonlinear scalar fields that are coupled 
to multiplets of vector and axial vector fields, and higher tensor fields as well. They are 
therefore difficult to test for everywhere regular solutions since their field structures, al- 
though mitigated by large symmetry groups, are very complex. Nevertheless it may be 
interesting to examine model theories obtained by coupling the dilaton field and a nonlin- 
ear scalar multiplet to the gravitational and electromagnetic fields, both for their intrinsic 
interest and also for their similarity to more fundamental theories. We shall consider a 
model theory with the following action which exhibits some of these features. 

S = J d 4 x^(R - KL) (1.1) 

N N (1.2) 

i=i i=i 

where the neutral scalar ip is coupled to as a dilaton, where e > 0, and where 
the charged scalars are coupled in the standard way to both the electromagnetic and 
gravitational fields by the following covariant derivative 

V ll = d l > + —A li + T li . (1.3) 

Here Y ^ is the gravitational part of the covariant derivative and V is the nonlinear scalar 
interaction. It would be possible to include magnetic couplings in (1.2) but we shall 
consider only the case in which magnetic field and magnetic charge vanish. We have not 
included any self-coupling of the dilaton field. 

When ifj Q is required to be constant, we have a theory with standard couplings. 

One may think of this model theory as interpolating between standard physics 
ipo = constant, and new physics associated with a Kaluza-Klein theory. It is not, however, 
intended as a truncation of a supergravity theory. 
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II. THE SOLITON PROBLEM 



In the action (1.1) there are the following coupled fields: the metric g^x, the vector 
potential A^, and the scalar fields tpi. The equations of motion are 



We require everywhere regular solutions of these equations subject to solitonic bound- 
ary conditions. Therefore no central singularities or horizons are permitted. These soli- 
tonic solutions concentrate energy- momentum and charge- current within lumps of field. 
In general, the mass spectrum of these particle-like lumps may be discrete. 

In this theory both the nonlinear scalar fields and the Maxwell field contribute to the 
source of the gravitational field while the complex scalars are the source of the electro- 
magnetic field. To avoid irrelevant complications, it is assumed that the charged scalar 
fields are spherically symmetric and harmonically time-dependent, while the dilaton is also 
spherically symmetric but time-independent: 



We shall need to consider only one member of the SU(N) scalar multiplet. 

Since we are dropping the magnetic fields, the vector potential also vanishes and the 
scalar potential is spherically symmetric. Thus the electric and gravitational fields should 
resemble the Reissner- Nordstrom solution except near the origin where the mass and charge 
are spread out and all fields remain finite. 




i = 0,...,iV 



(2.1) 



(2.2) 



(2.3) 



i/j (x) = Ro(r) 

ijj i (x) = e iu}t R i (r) , i = l,...N 



(2.4) 
(2.5) 
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III. THE FIELD EQUATIONS 



The field equations implied by (1.2) and (1.3) may be expressed as follows: 
(a) The gravitational equations 



R li x = KG liX K 
where k is Newton's constant. Here 



8ivk 



- ^ 1 T 

(b) The scalar equations 

* dV 

^ A v M v A ^ + ^ = o 

The operator V M contains the term in (3.5) but not in (3.4). 

(c) The electromagnetic equations 

The usual equation for the conserved current follows from (3.6), namely 

N 



To simplify the general equations (3.1), choose the line element 

ds 2 = e^c 2 dt 2 - e x ^dr 2 - r 2 (d6 2 + sin 2 6 dip 2 ) 
and in (3.4) and (3.5) choose 



F a p = 





—E 
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E 























Vo 








0/ 



where 



dp 



and ip is A , the scalar potential of the electromagnetic field. The complete set of gravita- 
tional equations is now 



v" 

T 


v'X 
+ 4 


(^) 2 
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v" 

T 


u'X' 
4 




) 2 


W 




= K 


y 



7/ 



= K 



e , ,. 

7 + ^l fl 



2 ^y(B 2 ) + a 



r 



±(R> ) 2 + (Rr + \e x V(R 2 )-A 



+ e A - 1 



The off-diagonal equations are satisfied identically. 
Here 

A = e eR °- v {<p') 2 /2c< 

N 

r 2 = J2 r "- 



(3.11) 
(3.12) 
(3.13) 



(3.14) 
(3.15) 



The equations for the dilaton and the charged scalars are: 





R o + - R 'o+ 2 



-R7 H — -R^ H — R'i + 



uj e 



+ — ip) e x -»R l -e x -—V(R 2 ) = 1 i = 1 . . . A(3.17) 
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Finally the single equation for the electromagnetic scalar potential is 



<p" + V 
r 



-(»' + \') + eR' 



2 nc\c nc 



The boundary conditions may be chosen as follows: 

^(0) = vA0)=0 
Ri(oo) = A(oo) = z/(oo) = y?(oo) = 
R'i(oo) = A'(oo) = z/(oo) = <f(oo) = 

where % = 0, . . . , N. In addition we require 

V(oo) = . 



(3.16) 



(3.18) 



(3.19) 
(3.20) 
(3.21) 



(3.22) 



Since the Ri are scalars, the equations (3.19) are invariant conditions which exclude 
cusps at the origin in these fields. 

The corresponding restriction on ip is 



ip'(0) = E(0) = 
5 



(3.23) 



which is required since the origin is at the center of a spherically symmetric charge distri- 
bution. In invariant form (3.23) reads 

(V A ) r=() = 0. (3.23)' 

In the rest of this paper we shall use the gauge where 

R i = 0, i = 2,...,N (3.24) 
Ri = R . 

We shall also need the following relation 
(v' + X) 



K 



(3.25) 



obtained by adding (3.11) and (3.12). It follows from the preceding equation and the 
boundary conditions that 

z/(0) + V(0) = . (3.26) 

The invariant form of (3.26) may be expressed in terms of g, the determinant of the metric 
tensor, as follows: 

U ] (3.27) 



9 flat / r —Q 

One also notes that (3.13) implies 

u'(0)-X'(0)=(-(e X -l)) ■ (3.28) 

Here 

\ = £ +£ 2 r 2 + ... (3.29) 
since A(r) and z/(r) are even by (3.11)-(3.13). If 

e A(0) ^ 1 (3.30) 

then by (3.28) 

i/(0) - A'(0) = oo (3.31) 

and by (3.26) both u(r) and A(r) would have infinite cusps at r = 0. If A(0) vanishes, 
however, then by (3.28) 

z/(0) - A'(0) = . (3.32) 
Therefore we adopt the following boundary conditions 

A(0) = (3.33) 
A'(0) = u'(0) = . (3.34) 
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These boundary conditions are compatible with (3.26) and (3.32) as well as with (3.11), 
(3.12) and (3.13). Finally by (3.25) we have 



p oo 

[u + A]!° = -K drr 
Jo 



1 



uj e 



R 2 



and since the integrand is positive, while K is negative, 

[* + A] ~>0. 

Then by (3.20) 



By (3.33) 



i/(0) + A(0) < 
i/(0) < . 



(3.35) 



(3.36) 

(3.37) 
(3.38) 
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IV. APPROXIMATE SOLUTIONS 



The boundary conditions (3.19) and (3.20) imply 

lim A(r) = (4.1) 

i — >o 

Therefore (3.11)-(3.13) all approach the corresponding equations (except for the frequency 
shift lo + e<fi/h) that hold in the absence of a charged source. These equations were 
previously studied. 4 Their solutions remain finite and have flat tangents at r = 0. 

(a) The Electromagnetic Scalar Potential. 

Let us consider (3.18) near r = 0. Guided by earlier results, 4 let us set 

R = a + a 2 r 2 + ... (4.2) 

R = a + a 2 r 2 + ... (4.3) 

A = £ 2 r 2 + £ 4 r 4 + . . . (4.4) 

v = no + n 2 r 2 + n^r A + . . . (4-5) 

Then 

R' (0) = R'(0) = Y(0) = t/(0) = (4.6) 
and (3.18) near r = becomes 



Set 



Now (4.7) becomes 



The solution of (4.9) is 



where 



Therefore 



v = - + f- («) 

c nc 



y" + 2 -y' = -\^~ e "° a ly- ( 4 - 9 ) 



sin pr , A _ _ 

y = (4.10) 



p 2 = -^e'^a 2 . (4.11) 



lim Lp = —= exp 



a 



ao--u . (4.12) 



Hence all fields including <p(r) remain finite with vanishing derivatives at r = 0. On the 
other hand at very large distances (3.18) becomes 



2 

<p" + -ip' = (4.13) 
r 

and 

ip = 9- (4.14) 

r 

where Q is the charge, 
(b) The Dilaton. 

Let us next consider the equation of the dilaton (3.16). Since all first derivatives vanish 
at the origin, it follows from this equation that Rq and R'q also vanish at this point. The 
first non-vanishing derivative at the origin is 

B™(0) = -2ee-^°) +ei?0 (°V'(0)) 2 • (4.15a) 

Let us first assume that e > 0. Then Ro(r) has a very flat maximum at the origin, and 
therefore near the origin 

R' (r) < (4.156) 
Equation (3.16) also implies that Ro(r) has no minimum anywhere since 

R'i(r) = -ee-"^^') 2 < (4.16) 

wherever R'o(r) vanishes. 

To determine the behavior of Ro(r) at very large r one may anticipate the result that 
A and v both decrease as 1/r. Then (3.16) becomes approximately 

^(r) + ^(r)+e(^) 2 = (4.17) 



and if one sets 



one finds by (4.17) 



R = Ar s + B (4.18) 

<j) = 9 (4.19) 
r 



Ro(r) = -^f + B (4.20a) 



To satisfy the boundary condition at infinity set B = 0. By (4.20a) one has at large r 

R' Q (r) = > (4.206) 
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Now (4.15b) and (4.20b) are in conflict unless there is an intervening minimum. But (4.16) 
does not allow a minimum anywhere, and so it is not possible to satisfy the boundary 
conditions at both r = and r = oo. If we assume the opposite sign of e, then there will 
be a minimum instead of a maximum at the origin and all the inequalities will be reversed; 
but there will remain the same incompatibility between the boundary conditions at the 
origin and at infinity. Since the presence of the dilaton field appearing in (1.2) therefore 
precludes the existence of a non-topological soliton, we shall set ipo = in (1.2) and delete 
Rq from all the equations following this paragraph. The resulting action and equations are 
then standard, but the potential remains general. 



(c) The Gravitational Equations. 

Let us next examine the set (3.11)-(3.13) at r = oo subject to the boundary conditions 
(3.20)-(3.21). (We shall anticipate the result that R vanishes faster than v and A as r — > oo.) 
Then by (3.25) we have at large r 

v' + X' = . (4.21) 

By (3.20) 

Therefore at large distances 
Then by (3.13) 



u(oo) + A(oo) = . 
u+X=0 . 



X' 



K 



+ 



e x - 1 



(4.22) 
(4.23) 

(4.24) 



where we have also substituted (3.14) and (4.23). 

By the boundary conditions V(R 2 ) vanishes at infinity. We know that A vanishes as 
1/r and it will be seen that R vanishes as a spherical Bessel function. Then (4.24) becomes 



A' K 

r ~ 2c 2 



(4.25) 



To second order in -. let us set 

T ' 



2m £ 2 
X = — + -§ 

~ Q , ?2 



(4.26) 
(4.27) 



Then by examining terms of order ^ one sees that 

KQ 2 



£ 2 - 2m z + 



2c 2 



Note that 



e , = (e _ A) = 1 _2^ + 2m^ 



(4.28) 



(4.29) 
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Then by (4.28) 

p" 1 2m (AIM 
€ =1 - — ~2^> (4 - 30) 

where to and Q are mass and charge. The terms in £ and agree with the Reissner- 
Nordstrom result for a point source. For the present case of an extended source, there will 
be additional terms of higher order in K 

(d) Charged Scalars. 

In order to extract similar asymptotic relations from (3.17), it is sufficient to assume 

V{R 2 ) = ^R 2 + ... . (4.31) 

Then at large r one has 

R" + ! R > + V -Z2L R > + (HL + w\ \^R - e x BR = . (4.32) 
r 2 \c nc J 

Let us carry R to the same order as A in -. Therefore, let us make the ansatz: 

R = —(l + P±) . (4.33) 



R" + ^R'= (a 2 + ^)R. (4.34) 



Then 

R " + r 

Eq. (4.32) may now be written to order ^R as follows: 

9 2p 1 a\ 2ma ( 4m 2£ 2 + 8m 2 \ uj 2 2e _ / 1 4m 

| 2eujg 2 | e 2 Q 2 / | 2m | l 2 + 2m' Q 
foc 2 r 2 h 2 c 2 r 2 \ r r 2 

From (4.35) one finds 



2 



C 

,2 



where 



By (4.36) and (4.37) the tail of R(r) is determined by 



« 2 = ^ + T-^-Q • (4-40) 



uo 2 iu 



c 2 hc 2 m 



If the mass is very large 



UJ 



(4.35) 



B = a 2 + -r > (4.36) 



2a; 2 ; e 
c 2 foe 2 m 

(pi - m)a + 2m 2 a 2 + #(Q) = (4.38) 



C 2 TlC 2 TO ^ ^ 



a = - . (4.41) 
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V. THE SPECIAL RELATIVISTIC EIGENVALUE PROBLEM 



The general relativistic equations that must be satisfied are (3.11)-(3.18) subject to the 
boundary conditions (3. 19)- (3.22) as well as (3.33) and (3.34), after the dilaton equations 
have been deleted. 

The qualitative nature of the solutions to these equations is already determined by 
the eigenvalue problem generated by the special relativistic limit of (3.17), namely 

r c z OR 

which holds in the limit at large distances. 

To discuss (5.1) introduce the auxiliary function: 

H= 1 -(R'f-V(R*)+ 1 -^ r R>. (5.2) 

By (5.1) 

dr r 

It follows that the representative point in the phase plane (the RLR'-plane) will always 
move towards lower values of H unless R' = 0. Hence, if H (R, R') is pictured as a surface 
above the -R-R'-plane, the representative point will fall toward minima on this surface. 
The extrema in the Redirection are given by 

This locus is the i?-axis. It is a minimum in the R' direction since 

92R =1>0. (5.5) 



A2 



d(R') 

The extrema in the Redirection depend on V(R 2 ). 

For definiteness assume that V(R 2 ) may be represented by the polynomial 

V(R 2 ) = |i? 2 + ^ 4 + ^ (5.6) 

where 

S>0, C<0, D>0 . (5.7) 
Since D is positive, the field energy is bounded below. Then 

i)H % + ^tR = —(BR + CR 3 + DR 5 ) + ^R (5.8) 



OR OR c 2 v ' c 2 

g| = -g + ^=-( B + 3^ + 5^)4. (5.9) 
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By (5.8) the extrema in the indirection lie at the roots of 

dH 



dR 

where 

The five roots of (5.10) are 
and 



— (B + CR + DR)R = (5.10) 



B = B-%. (5.11) 



R = (5.12) 



C\+ E' 



1/2 




Here 



Then 



R-+ = -R ++ (5.13) 

= -R+- (5.14) 

E 2 = C 2 — ADB . (5.15) 

-(SL-- 



We now require B = B- u 2 /c 2 > 0. Then the origin is a maximum in the .R-direction. 
At -R++ and _R_| we have 

fLhH " -l(\ C \ + E)<0 (5.17) 



V dR 2 J ++ D 

O + _-|(' c '- £ » >0 < 6 - i8 » 

since \C\ > E. 

Hence the points (R^ , 0) and (—R-\ , 0) are minima in both R and R' directions. 

The point (0, 0) is a minimum in the R' direction but a maximum in the i?-direction and 
is therefore a saddle point. Similarly the points (R++, 0) and (— R++, 0) are saddle points. 

The phase portrait of H(R, R') is shown in Fig. 1. We shall assume in this figure that 

H(R ++ ,0) > #(0,0) or H(R ++ ,0)>0 (5.19) 
To satisfy this condition, D is restricted by the following inequality: 

3 C 2 

- T >D. (5.20) 
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If V(R) were a higher order polynomial or some other functional, the phase portrait 
could have more structure than is shown in the figure. It will become clear, however, that 
only that portion of the surface near the origin is relevant here. 

Let Co and C++ be iso-H curves passing through (0,0) and (R++, 0) . Require that 
D be restricted by (5.20). Then if the initial point of the solution curve lies anywhere in 
the phase plane between Co and C++, this curve must terminate, according to (5.3), at 

(0,0), (.R-i , 0) or (— -R-i , 0). In the soliton problem the initial point of the solution curve 

must lie on the .R-axis, and if it also lies between Co and R++ it must terminate at either 
one of the two attractors, R-\ and —R-\ , or at the origin. 

In Fig. 2 the region between Co and C++ is expanded. The solution curves are shown 
starting at po, so, and qo and ending at P, S, and Q. In Fig. 3 the same solutions are 
shown in configuration space. 

Every solution curve starting on the axis between R^ and R++ and, depending on 

its initial point, lies in one of two classes going to either the left (P) or right (S) attractor, 
unless it goes to the origin. The relation between curves in the p and s classes and the 
eigensolution is 

~R(p) + R(s) 



lim 

Po—S - 



= R(q) (5-21) 



2 

where p, s, and q refer to the three curves going to P, S, and Q. 

By bringing p a and s a together one may determine the eigenvalue q Q with arbitrary 
accuracy. 

Similarly Eq. (5.21) is a prescription for obtaining the eigensolutions R(q) with arbi- 
trary accuracy. 

Note that at large r we have approximately 

R(p) = A + u(p) (5.22) 
R(s) = -A + u(s) (5.23) 
R(q) = u(q) (5.24) 

where A = R + -. Since (±A, 0) are constant solutions of (5.1), they are also roots of the 
following polynomial and hence extrema of H: 

Here u p , u s , and u q are spherical Bessel functions that solve (5.1) after it is linearized near 
±A and zero. The linearized equation is 

U " + l U ' + (Df) U = ° ^=(±A0) (5.26) 

and 

ci 

u(p) = — sin(cl> p r + e p ) (5-27) 

u(s) = — sin(u) s r + e s ) (5.28) 

u(q) = %" ar (5.29) 
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where 




d 2 H 



) 



(5.30) 



OR 2 



A 



) 







(5.31) 



The functions u(p) and u(s) describe the spiral motion toward the attractors while 
u(q) describes the exponential approach to the saddle point. Note that the oscillations of 
u(p) and u(s) begin to cancel in the sum ^(R(p) + R(s)) because of the phase difference 
between the arguments of u(p) and u(s) that arise because the motions around S and P 
are clockwise and counterclockwise respectively (u p = —uj s ). The phase difference for the 
actual numerical solutions is shown in Fig. 4a. In the large distance limit the oscillatory 
behavior disappears since \{u{p) +u(s)) must approach u(q). 

As p and s Q are brought closer to q a the oscillatory behavior of p and s recedes to 
infinity and both p and s approach the eigensolution, q, which terminates at the origin 
with no oscillations. In Fig. (4b) this eigensolution is shown in configuration space. 

The preceding discussion continues to hold as the initial values, p and s, are increased. 
It may then happen that the p and s curves cross the R' axis several times before being 
caught on one of the attractors. In that case the eigensolutions will node a corresponding 
number of times and the resultant field structure may be described as an excited state of 
the soliton. 

On the other hand, if the initial point lies inside the separatrix, the solution curve will 
always end at the nearby attractor. According to one of the boundary conditions that we 
have imposed, namely R(oo) = 0, such a curve is not an eigensolution. If the boundary 
conditions were relaxed, however, so that R(oo) is interpreted as a vacuum expectation 
value, as in the case with gauge solitons, this interior curve could be regarded as belonging 
to a continuum of eigenstates. 

The preceding qualitative discussion is based entirely on equation (5.1) but the illus- 
trating figures are based on the numerical solution of the complete set of simultaneous 
differential equations. This is possible because the qualitative behavior of the solution to 
(5.1) is very similar to the corresponding solutions of (3.17) viewed as one member of the 
set of simultaneous differential equations. 
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VI. THE GENERAL RELATIVISTIC EIGENVALUE PROBLEM 



Although we have not found an if-functional appropriate to the full set of equations, 
nevertheless the numerically obtained solutions of the original equation (3.17) behave qual- 
itatively in the way just described for the simplified equation (5.1). 

We may gain additional information by regarding the full equation (3.17) as a pertur- 
bation on (5.1) and by confining the discussion to the RR'-p\a,ne. 

Consider, for example, the region of phase space where v = X = v' = X' ~ 0. Then 



R" + -R' + 

r 



and 



dH 

dr 



= R' 

= {R'f 



u e 

c hc^ 



dV 



R- 



dV 
OR 







(6.1) 



R" + u z R - — 



2 

r 



OR 

c Tic^ 



R_ 
R 1 



(6.2) 



For very large R 



if 



dH 

dr 



2(R 



t\2 



-1 + 



ue Q 
he 2 a 



Q 



and R 



Then 



dH 
hm — - — = 

r^oo dr 



(6.3) 



(6.4) 



(6.5) 



but it may be of either sign depending on the sign and magnitude of the second term. 
If the sign becomes positive after the representative point crosses the separatrix, it will 
settle into a closed orbit about either the left or right minimum, depending on whether 
it is caught in the left or right lobe. Then in configuration space R will not approach a 
constant at infinity but will continue to oscillate. 

We also note that the sign of dH/ dr may be switched outside the separatrix by slightly 
changing the magnitude of <p. If this happens when the representative point has come close 
in the approach to either lobe, it may happen that it finally terminates in the opposite 
lobe. In configuration space, this behavior will appear as a discontinuous change in the 
solution curve R(r) caused by a small change in the parameter (p(0). 

None of these remarks invalidate the exact, non-perturbative procedure for determin- 
ing the eigenvalues and eigensolutions of the set (3.11)-(3.18) (always excluding the dilaton 
equation). The most important result of the numerical work is the observation that the 
R(r) solutions divide into two classes when the four general relativistic differential equa- 
tions are simultaneously integrated. This key feature of the special relativistic problem is 
exactly preserved. Now, however, associated with R(p), R(s) and R(q) are the families 
[X(p), v{p)-> <p(p)], [A(s), f(s), <p(s)] and [X(q), v(q), <p(q)}- Then in addition to (5.21), 
we have 
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lim 

Po— so— *-o 



A(p) + A(s) 



= A( ? ) (6.6) 



with similar relations for z/ and <p, so that in the limit [R(q), A(g), z/(<z), <£>(<z)] is a complete 
set of solutions corresponding to q . We also note that the p and s curves for A, v and ip 
are nearly the same while the corresponding curves for R are of course quite distinct at 
large r. 

There is also the following approximate procedure: Carry the integration into the 
region where the gravitational potentials become asymptotic. Up to this point, say f, 
which is somewhat arbitrary, approximate the four functions by their s and p average: 
\[R{p) + R(s)], |[A(p) + A(s)] — We shall choose as a criterion for determining f the 
condition: 

A(f) + u(r) = . (6.7) 

Beyond this point approximate the eigensolution R(q) by its asymptotic form (a spherical 
Bessel function). Then match the asymptotic form of R(q) to ^[R(p)+R(s)] at f as follows: 

L (R(p)+R(s)) f (6.8) 



f 2 

{^^)=\{R'{P) + R'{s)) f . (6.9) 

To complete the integration beyond f feed the so obtained R(q) into the complete set of 
differential equations while taking ^ [A(p) + A(s)] ... as initial values of the other functions. 
If the point defined by (6.7) is difficult to reach, one may instead choose f as the point at 
which the curves s and p begin to separate. Figure 4b was obtained by this alternative 
procedure. 

To find the functions (f(r), A(r) and z/(r), corresponding to the eigenfunction R(r), 
we have integrated backwards by feeding R(r) into the appropriate differential equations 
while imposing initial conditions determined by the asymptotic forms (4.30) and (4.31) 
and computed for a point and a mass m. The values of and m must be adjusted 
to agree with the boundary conditions at r = 0. The results of this backward integration 
are presented in Fig. 5. 

Additional checks may be obtained by integrating equations (3.13) and (3.25) as fol- 
lows: 

-j poo poo 

[u(r) - A(r)]~ = - / dr(e x -l)+K drr[V(R 2 )e x + 2A] = J x (6.10) 



poo 

[v(r) + A(r)]!° = ~K drr 

Jo 



\ 2 to2 



= J 2 ■ (6.11) 



Now imposing the boundary conditions v{po) = A(oo) = we have 

A(0) - i/(0) = Ji (6.12) 
A(0) + i/(0) = - J 2 < (6.13) 
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These equations may also be used iteratively to improve approximate solutions already 
obtained. 

The discrete spectrum of .R-eigensolutions will induce discrete spectra on the other 
fields. As a consequence the mass of the soliton, which is identified by the coefficient of the 
1/r term in the asymptotic gravitational potential, is also limited to a discrete spectrum. 

It has been noted that a small change in </?(0), like a small change in R(0), may cause 
the solution curve in the .R-R'-plane to jump between the P and S lobes. Therefore we may 
divide the [<£>(0), -R(O)] plane into + and - areas according to whether the solution curve 
terminates in the P or S lobe. More generally, one may divide the 4-dimensional space 
of initial values with coordinates [A(0), z^(0), (p(0) 7 -R(O)] into + and - subspaces. Then 
any boundary between a + and - subspace will describe a hypersurface of eigenvalues. We 
have investigated only a limited region of the 4-dimensional parameter space. 
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VII. DISCUSSION 



General relativistic solitons have also appeared in the literature on boson stars. 5 The 
purpose and methods of this paper are quite different, however, since they are more di- 
rectly related to the first reference 1 and to the possible existence of non-topological soli- 
tons in fundamental completions of Einstein theory. The main result of the present work 
is the demonstration that the method of analysis employed in the earlier work may be ex- 
tended to many coupled fields. While the field structures studied in this paper couple the 
gravitational, electromagnetic, nonlinear scalar, and dilaton fields, there are many other 
interesting possibilities that are tractable by the same method, such as gauge solitons and 
other structures generated by Kaluza-Klein-like theories. 

Since there has been a particular emphasis in our investigation on the nonlinear scalar 
potential, the generality of this potential is of interest. Although a particular polynomial 
representation was chosen, it is clear that the analysis depends only on the distribution 
of minima and saddle points of the dissipative function along the field-axis in phase space 
and therefore on the distribution of maxima and minima in the potential. One also notes 
that there are essentially only two kinds of potential, namely Higgs and inverted Higgs, or 
volcano, which are related by the interchange of minima with saddle points. If the potential 
is of the volcano type, there is a saddle point in the dissipative function flanked by two 
minima and if it is of the Higgs form, there is a central minimum bracketed by two saddle 
points. If this pattern of three extrema is centered at r = 0, then the eigenstates of the 
former vanish at great distances. In the Higgs case there is a continuum of eigensolutions 
and in the case of the volcano potential there is a discrete set of eigensolutions, as we have 
discussed. If these patterns are not centered at R = 0, however, then the corresponding 
eigenfunctions do not vanish at infinity. In this case one has the situation encountered in 
gauge theories where the limiting value at infinity is interpreted as a vacuum expectation 
value. If the potential is general with several minima, there should be both continuous 
and discrete spectra associated with it. For example, the potential may be Higgs-like at 
weak fields and volcano-like at strong fields. If the boundary conditions are interpreted 
in terms of vacuum expectation values, then the physical constants appearing in the field 
Lagrangian are encoded as vacuum expectation values. 

Since our remarks concerning the dilaton field may also be relevant to certain string 
models, let us consider two examples. ( 6 ' 7 - ) 

A dimensionally reduced superstring in four dimensions can be described in terms of 
N=4 supergravity: the action associated with the SO(4) version of this theory is 

-R + 28^8^ - (e'^F^F^ + e 2 ^G^G^)] (7-1) 

as expressed for example in reference (6). Dropping the G term one sees that the dilaton 
field (j) appears here in the same way as in our equation (1.2). Therefore we may conclude 
that the same conclusion also holds here, namely, that non-topological solitons do not exist 
for this theory. 

As a second example consider the 5-brane soliton that is a solution of the field equa- 
tions describing the low energy heterotic string. ( 7 ) This soliton embodies a 4-dimensional 



I{SO{A)) = / d 4 x 
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Yang-Mills instanton coupled to dilaton and axion fields according to superstring theory. 
The dilaton equation in this case is^ 7 ) 



= ±—a'e^TrF^F p<J (7.2) 

This equation also resembles our equation (3.4) but with the following differences: the 
source term in (7.2) is pseudoscalar rather than scalar and moreover it lies in a gauge alge- 
bra; the instanton exists not in Minkowski but in Euclidean space, and the dilaton at spatial 
infinity does not vanish but approaches its vacuum expectation value. The field of this 
topological soliton is everywhere finite as required and its extension is determined by the 
size of the instanton. The full structure has not been presented, however, and its interpre- 
tation in conventional spacetime depends on how the 5-branes wrap around topologically 
non-trivial surfaces in the internal compactification manifold. Since this topological struc- 
ture is is therefore essentially different from the non-topological four dimensional solution 
studied here, our remarks about the dilaton field do not apply to this example although 
they do to the preceding case. On the other hand, the procedure described here may be 
useful for investigating non-abelian as well as abelian theories, and therefore for describing 
topological as well as non-topological solitons. 
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FIGURE CAPTIONS 



Fig. 1. Phase portrait of H(R, R'). Q, Q' and Q" are saddle points. S and P are minima. 

The R coordinates of Q, Q' and Q" are R |_, and R++ respectively. The R coordinates 

of S and P are R^ and R respectively. 

Fig. 2. Solution curves s and p bounding the eigensolution q, shown in the phase plane. 

Fig. 3. Solution curves s and p bounding the eigensolution q, shown in configuration space. 

Fig. 4. (a) The asymptotic behavior of the solutions s and p, illustrating the approximate 
phase cancellation at large r. (b) An eigensolution of R(r) in configuration space as 
determined by the numerical procedure described in the text. 

Fig. 5. (a) The gravitational eigensolutions A(r) and v{r). (b) The electromagnetic 
eigensolution <p(r). 
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